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I.  INTRODUCTION 


The  trajectory  equations  that  govern  the  kinematic  behavior  of 
missiles  flying  various  proportional  navigation  trajectories  are  pre¬ 
sented  with  the  goal  of  establishing  missile  maneuvering  requirements 
for  a  specific  set  of  engagement  parameters.  These  parameters  are 
those  that  might  be  encountered  in  ground-to-ground  tank  engagements 
at  the  2  to  3  kilometer  range.  Analog  simulation  results,  based  on  the 
trajectory  equations,  show  the  maneuver  capability  required  of  a  homing 
missile  to  overcome  various  perturbations. 

In  addition,  a  problem  peculiar  to  seekers  that  track  the  target 
image  contour  is  treated.  This  problem,  that  of  guidance  loss  in  the 
terminal  tracking  phase  due  to  the  increase  in  image  size  with  decreas¬ 
ing  range,  results  in  a  need  for  the  missile  to  coast  unguided  for  a 
fraction  of  this  trajectory.  The  theoretical  miss  distances  due  to  this 
phenomenon  have  been  calculated  and  are  also  presented. 

II.  MISSILE  PERFORMANCE  REQUIREMENTS  FOR  PROPORTIONAL  NAVIGATION 

In  proportional  navigation,  the  missile  turning  rate  (  y  )  is  made 
proportional  to  the  line-of-sight  (LOS)  rate  (a)t 

y  ■  No- 

N  is  the  proportionality  (navigation)  constant.  Figure  I  shows  that 
the  normal  velocity  component  of  the  missile  is  given  by:  y  *•  V  ainy 
» Vyfor  small  angles;  therefore,  the  normal  acceleration  1b  given  by: 
y  n'Vy  +  Vy  »Vy  for  a  constant  velocity  missile. 

Thus: 

y  B  and  y  B  NVo- 

One  practical  means  of  attaining  proportional  navigation  is  to 
carry  onboard  the  missile  a  tracking  device  capable  of  measuring  the 
LOS  rate  and  using  this  information  to  command  a  lateral  acceleration 
proportional  to  this  rate.  A  means  of  measuring  <r  can  be  explained 
with  the  aid  of  Figure  I.  A  gimballed  seeker  attempts  to  track- the 
target,  i.  e.  reduce  the  error  angle  *  to  zero.  This  angle  exists  be¬ 
cause  the  tracker  cannot  have  an  infinitely  fast  response  and  will  thus 
lag  behind  the  LOS  by  an  amount  determined  by  the  tracking  loop  time 
constant  (T  ).  Thus  the  instantaneous  value  of  the  angle  p.  is  given 
by: 

^  =  J+TS 


From  Figure  I  it  is  obvious  that*  a<r-ji ,  and  thus: 
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-  -o-  _  — g_  .  <r  4-g-TS-o-  .  .js__  g 
I  +vS  1+rS  I  +  tS 

Therefore,  the  instantaneous  tracker  error  angle,  a  rather  easily 
measured  quantity,  can  he  used  to  effect  proportional  navigation  if  we 
multiply  it  by  some  gain  factor  X  and  use  the  resulting  value  to  com¬ 
mand  a  missile  acceleration: 


a"> ' *  T+TT  '  <» 


where  a  «=  missile  lateral  acceleration  in  g's, 
m 

The  above  equation  will  be  used  in  the  derivation  of  the  homing  tra¬ 
jectory  differential  equations.  These  equations  are  greatly  simplified 
if  we  use  the  constant-bearing  (collision)  course  (Reference  1  and  2) 
as  a  reference  and  consider  the  proportional  course  as  perturbations 
about  this  reference.  Such  a  reference  frame  is  shown  in  Figure  IX, 
and  the  associated  symbols  are  identified  in  Table  I.  From  the  figure, 
we  obtain: 

Zm  £  ~am9  (2) 

i.  e.  the  commanded  acceleration  of  the  missile  normal  to  the  reference 
course  is  directed  so  that  it  tends  to  eliminate  Zm,  thus  the  negative 
sign.  Combining  equations  1  and  2: 


g\t 

1+rS 


or,  converting  the  whole  equation  to  time  domain: 


r 
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•••  « • 

rZm  +  Zm  *  -gXrd- 


Integrating  above : 

TZm  +  Zm  =  -fl\TO-+C 


(3) 


C  is  the  constant  of  integration  (initial'  conditions)  and  is  given  by: 

c=TZm|0*  +  Zm|0+  +gXr(o-|0J 

where  o+  denotes  value  of  variable  at  t  a  o 
Equation  3  can  then  be  written: 
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(4) 


vZm  +  Zjns-gXrcr  4-  rZm|o+  +  2m|o+  +gXr(tr|04) 

Equation  k  is  then  the  trajectory  differential  equation  and  can  be 
solved  if  <r  can  be  determined.  Sincecr  is  generally  small,  the  sim¬ 
plifying  approximation  cr  «  tan  «r  can  be  used.  Prom  Figure  II  it  is 
also  clear  that  the  tangent  of  a  is  given  by: 

tnn  rj  -  (T+2M3+4) 

R  t 


(Xmsin8  -f  ZmCosS)  -  (Xt  sin<M-Ztcosfr) 


Since  R(.  *»  R  (tQ  -  t)  where: 

«  a  range  rate  ■  Ve  ■  closing  velocity 
t  n  nominal  homingctime  « 
t°  a  elapsed  flight  time 

Introducing  above  value  of  a  into  equation  k  and  substituting 
Vc  (t  -  t)  for  R,  we  obtain  the  complete  trajectory  differential 
equation : 

TZm  +  Zm«  y^^(Xmsin8+Zmco88-Xtsim//-Ztcos^)+TZm|0++Zll,Jtt4 

+gXr(o-|ot) 

Transposing  cos  8  to  the  left  side,  we  obtain: 

T2n+Zm+^^°s®j  Zm=  }  (Xmsin8-X, sin^-ZfCOSVr)  (5) 

9XtCos8  +  rMo*  +Zm|0+  J 

The  term  yc  is  almost  constant  for  small  deviations  from  a 

collision  course  and  is  generally  defined  as  the  effective  navigation 
constant  N'.  A  more  common  expression  for  N'  is  given  by: 

m'b  Vm  Nc°s8 

Vc 
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The  above  two  expressions  can  be  shown  to  be  equivalent  if  the  follow¬ 
ing  equality  can  be  demonstrated : 

flXr  *VmN 

This  can  be  done  as  follows: 

It  has  been  shown  that  the  seeker  measures  r&  but  with  a  lag 
(I+ts)  term.  Thus,  since  the  basic  error  is  rcr  ,  then  from  equation 
1,  the  basic  commanded  acceleration  is: 

dm  "  Xro-  (6) 


From  the  basic  navigation  law,  we  have  shown: 

yaNVmo-  or  ama— 2-d-  (7) 

Equating  the  right  sides  of  equations  6  and  7,  we  obtain: 

Xr°NVm/g  or  Xrg*NVm 

Dividing  and  multiplying  the  right  side  of  equation  5  by  cos  8  and 
substituting  N',  we  obtain: 


t  Zm  -fZfn  + 


N*  Zm 

♦o“  t 


— (Xmtan8-Xt 

To  ' 


simfr  cos^ 
cos  8  1  cos  8 


+  ‘rZ|o++gXr(<r|o^)i-2|ot 


(8) 


The  preceding  trajectory  equation  is  a  linear  differential  equa¬ 
tion  (time  varying  coefficients).  Since  the  principle  of  superposi¬ 
tion  applies,  the  equation  can  be  solved  for  certain  specific  condi¬ 
tions  of  interest  one  at  a  time  and  the  results  linearly  combined.  Of 
particular  interest  is  the  solution  of  Zm  for  the  following  condi¬ 
tions  : 

1.  An  initial  missile  heading  error,  i,  e,  the  missile  is  not 
pointed  at  the  target  at  launch. 

2.  An  acceleration  bias.  This  can  be  caused  by  thrust  malalign¬ 
ment,  improper  gravity  correction,  etc. 

3.  A  maneuvering  target. 
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If 


An  initial  tracking  error,  such  as  would  be  caused  by  improper 
seeker  pointing  at  launch. 


Given  typical  values  for  the  above  four  conditions,  the  solution  of 
equation  8  for  2  m  yields  the  important  specification  of  required 
missile  maneuverability,  The  analog  computer  studies  to  be  discussed 
are  concerned  with  obtaining  this  specification,  since  no  one  specific 
missile  is  under  consideration,  a  range  of  representative  parameters 
have  been  selected  so  that  the  data  may  be  applicable  to  a  wide  range 
of  vehicles.  The  general  environment  considered  is  that  of  a  ground-to- 
ground  engagement  at  ranges  of  6000  and  10,000  feet.  Since  in  this  case 
gravity  acceleration  acts  normal  to  the  flight  path,  the  results  obtain¬ 
ed  are  valid  for  either  pitch  or  yaw.  This  is  due  to  the  fact  that  for 
the  ground-to-ground  case  an  external  "g"  bias  would  be  required  to 
overcome  the  gravitational  acceleration,  and  is  thus  not  part  of  the 
closed  loop  dynamics.  If  the  navigation  law  were  to  be  used  to  over¬ 
come  the  effect  of  gravity,  an  excessive  portion  of  the  missile's  ma¬ 
neuvering  capability  would  have  to  be  employed;  additionally,  the 
ground  clearance  problem  would  become  overwhelming  unless  the  missile 
were  to  be  launched  into  a  high  initial  trajectory.  This  cou}.d  be  done 
by  initiating  guidance  after  launch,  and  the  problem  would  be  reduced 
to  one  of  overcoming  an  acceleration  bias,  a  case  treated  herein. 


A.  Case  1  -  An  initial  heading  error 

From  Figure  II  it  is  evident  that  Z  m  1 0  +  =  VmS  in  yo  »  Vmyn 
With  all  other  initial  conditions  »  0,  equation  8  becomes: 


z+zm  + 


N‘Zm 


to-t 


=  VmXc, 


or,  since  f  -t  «R/Ve 


7  4.7  4.  N'Zm 

rZm+Zm  +  y 


*“  Vmy0 


(9) 


A  computer  program  for  equation  9  is  shown  in  Figure  III.  The 
scale  factors  shown  ere  only  one  set  of  several  used.  The  scale  changes 
were  necessitated  by  the  dynamic  range  of  the  equations  and  by  noise 
problems.  Graphs  1  through  9  show  plots  of  2m  versus  range  for  a  set 
of  representative  conditions.  (^>  in  the  graphs  refer  a  to  y0  ). 

Several  interesting  observations  can  be  made  from  these  plots.  The  in¬ 
dividual  plots  on  Graphs  1  ancL  2  show  the  effect  of  changing  r  only. 

The  magnitude  of  the  required  2 rn,does  n0^  c^anSe  with  T  ;  however,  the 
range  at  which  the  peak  value  of  2m occurs  does  change  considerably. 

This  is  also  evident  from  Graphs  4  and  5  where  the  initial  conditions 
are  the  same  as  for  1  and  2,  except  for  initial  range.  The  longer 
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range  results  In  a  decrease  in  required  peak  acceleration. 

Graphs  3  and  6  through  9  show  the  effect  on  the  required  Z-,  of  in¬ 
creased  N'  and  velocity.  It  is  evident  that  an  increase  in  either  para¬ 
meter  results  in  an  increase  in  missile  maneuver  requirements.  The  im¬ 
plication  is  that  missile  flight  time  and  navigation  ratio  are  the 
critical  parameters  of  interest,  This  is  evident  from  equation  9.  If 
we  replace  Vc7R  by  l/ttf-t;  Jerger  (Reference  1)  has  shown  that  for  a 
given  navigation  constant  and  initial  conditions,  the  acceleration  re¬ 
quirements  are  dependent  solely  on  the  ratio t0/r  .  This  term  is  de¬ 
noted  "control  system  stiffness.  "  For  example,  in  the  first  set  of 
plots  in  graph  l,to/T  *  12.  The  2m  values  would  be  equally  valid 
for  R  o  12,000  feet  and  V  =  2000  ft/sec  since  \q/t  would  remain  unchang¬ 
ed.  Another  useful  generalization  can  be  stated  as  follows.  For 
tb/r  >^Q  ,  the  system  does  not  deviate  considerably  from  a  perfect 
system  (fo/r»flO  ). 


B.  Case  2  -  An  acceleration  bias 


If  an  acceleration  bias  exists  in  the  missile,  an  additional  fix¬ 
ed  acceleration  will  be  commanded  in  addition  to  that  called  for  by  the 
navigation  law.  Equation  1  then  becomes: 


•m 


Qb  + 


XtS 

1+  Ts 


IT 


The  term  will  appear  after  one  integration  as  a^t  in  equation  3  and 
as  -a^t  on  the  right  side  of  equation  8. 

If  we  assume  no  perturbations  and  no  initial  heading  error,  equa¬ 
tion  8  reduces  to: 


TZm+  zm  + 


Since  2  =  -  Q. 

equation  reduces  to:  " 


and(t0-t)=R/V0 


the  above 


rZm  +  Zm  + 


N'V 

R 


~  Zm,=  —  aJt+T) 


(10) 


The  right  side  (forcing  function)  of  the  above  equation  is  the  only  term 
that  differs  from  equation.  9.  Thus,  with  the  one  exception,  the  same 
computer  program  may  be  used.  Mechanization  of  the  term  -0m(f  +  T ) 
is  shown  in  Figure  IV.  Plots  of  the  values  of  Z  fn  required  to  overcome 


6 


bias  errors  of  1 /k,  l/2,  and  1  "g"  are  shown  in  Graphs  10-16. 

C.  Case  3  -  A  maneuvering  target 

The  case  for  a  maneuvering  target  can  be  considered  as  follows: 

^  -M,,  =H*  H„.=x»=x*  =° 

and  Zt  =  |attJ  (constant  acceleration  target). 

Then  equation  8  can  be  written: 


rZm  +  Zm  + 


m 


N'VC 

R 


‘  2m  = 


N'VC  COS  If/  Qt  tz 
2  R  cos  S 


If  we  consider  the  most  severe  case,  i.e.  the  Initial  missile- 
target  collision  course  is  along  a  straight  line  (head  on  engagement) 
and  the  direction  of  target  maneuver  is  normal  to  this,  the  above  equa¬ 
tion  can  be  written: 

7  +7  +  N'VcZm  =  _  N'Vc  at  t2 

r '  *-m'  d  “  op 


This  equation  is  again  identical  to  the  other  cases  considered  except 
for  the  forcing  function.  Thus  the  computer  program  on  Figure  III  can 
again  be  used  with  a  suitable  modification  for  the  term  of  the  right 
side.  This  is  shown  in  Figure  V.  Plots  of  the  required  missile  accele¬ 
ration  for  two  magnitudes  of  target  maneuver  (4  ft/sec2  and  8  ft/sec2) 
and  a  representative  set  of  engagement  parameters  are  shown  in  Graphs 
17  through  22.  From  these  graphs,  it  is  evident  that  the  missile  maneu¬ 
ver  requirements  increase  considerably  as  the  target  is  approached. 


D.  Case  4  -  An  initial  tracking  error 

If  the  missile  borne  seeker  exhibits  a  tracking  error  at 
t  =  0,  due  to  imperfect  aiming,  it  will  deviate  from  its  correct  course 
and  restoring  commands  are  generated.  The  effect  of  can  be  eval¬ 
uated  by  solving  equation  8  for  the  case  where  it  is  the  only  perturba¬ 
tion.  It  has  been  shown  (Reference  1)  that  the  trajectory  equation 
forcing  functiort  for  this  case  is  given  by:  —  VCN  <0/COS  S 

Thus  the  trajectory  equation  becomes: 


rZm+Zro  +N'VcZm/R  =  VsN'ce/cosS 


The  mechanization  of  the  forcing  function  is  simply^  properly  scaled 
step  input  and  is  shown  in  Figure  VI.  The  plots  ofZmfor  various  para¬ 
meters  of  interest  are  shown  in  Graphs  23  through  29. 
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The  plots  indicate  that  very  large  initial  accelerations  are 
required;  however,  we  should  keep  in  mind  that  the  values  of  «0  employed 
are  also  quite  large.  Initial  tracking  error  much  smaller  than  1°  -  3° 
should  be  realizable  without  excessive  difficulty. 

IIL  THEORETICAL  MINIMUM  MISS  DISTANCE  VERSUS  COAST  RANGE 

The  employment  of  an  image  contour  tracker  as  the  error  sensing 
element  in  a  homing  missile  presents  a  multitude  of  peculiar  problems. 
One  of  these  is  related  to  the  angle  subtended  by  the  target  in  the 
field-of-view  and  is  of  interest  due  to  the  following  reasons: 

1.  At  maximum  range,  the  field-of-view  of  the  optics  must  be  small 
enough  to  yield  a  target  image  large  enough  to  exceed  the 
resolution  limitations  of  the  TV  seeker. 

2.  As  the  missile  nears  the  target,  the  image  will  subtend  a  pro¬ 
portionally  larger  portion  .of  the  field-of-view,  until  even¬ 
tually  the  image  to  field  area  ratio  becomes  unity  and  track¬ 
ing  information  is  lost.  The  missile  must  then  coast  unguided 
to  the  target.  For  hardware  presently  under  development,  the 
field-of-view  limitation  is  approximately  2°  for  an  initial 
range  of  2  kilometers.  Considering  a  7  1/2  x  7  1/2  foot  tar¬ 
get,  the  missile  must  coast  for  approximately  the  last  250 
feet.  For  a  3  kilometer  initial  range,  the  coast  distance  be¬ 
comes  approximately  500  feet. 

Miss  distance  calculations  are  extremely  complex  and  involve  many 
probabilistic  and  design  factors.  Our  concern  here  is  to  determine  the 
misa  distance  due  only  to  target  maneuver  during  the  missile  coast 
period  when  trajectory  corrections  cannot  be  effected.  Two  alternatives 
are  discussed. 

A.  Neutral  Command  Terminal  Trajectory 

In  this  scheme,  when  loss  of  track  occurs,  the  control  sur¬ 
faces  return  to  neutral  and  result  in  a  straight  line  missile  trajectory. 
The  miss  distance  can  then  be  approximated  by  the  target  travel  normal 
to  the  collision  course  during  the  missile  coast  period.  These  values 
are  tabulated  in  Table  I  for  two  values  of  target  acceleration  and 
initial  range,  and  for  several  values  of  missile  velocity.  The  miss 
distance  (m)  is  given  by: 

m  o  1/2  attg 

where  at  =  target  acceleration  normal  to  LOS 
ts  =  coast  flight  time 

B.  Retention  of  Last  Command  Prior  to  Loss  of  Track 


If  the  control  system  is  locked  in  a  certain  position  at  some 
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time  t,  the  miss  distance  for  maneuvering  target  can  be  found  by  study¬ 
ing  an  equivalent  case.that  of  acceleration  saturation.  This  occurs 
whan  the  engagement  dynamics  due  to  target  maneuver  are  such  that  the 
navigation  law  calls  for  a  command  that  exceeds  the  missile's  capability. 
The  two  cases  are  analogour,  and  the  miss  distance  can  be  calculated  by 
assuming  that  acceleration  saturation  occurs  when  the  seeker  loses 
track.  The  derivation  of  m  is  given  in  Reference  1  for  the  case  where 
T*  0: 


m  *  (t0  -  ts)  N'  x  t02 

2 

where  t0  =  R/Vc 

tB  -  tQ  -  R/Vc 

Miss  distances  were  calculated  for  the  same  parameters  as  Case 
A  and  the  results  are  shown  in  Table  II.  It  is  evident  that  the  values 
of  m  are  negligible.  It  is  true  that  the  r  -  0  assumption  is  not  quite 
valid}  however,  an  useful  approximation  is  obtained.  Additionally,  one 
general  important  conclusion  can  be  made,  i.e.  it  is  at  least  theoreti¬ 
cally  possible  to  hit  a  moving  tank-sized  target  with  a  missile  that 
employs  the  contrast  controut  tracker  (CCT)  under  consideration 

IV.  CONCLUSIONS 

It  has  been  ohuwn  that  for  a  perfet  system,  the  miss  distances  re¬ 
sulting  from  terminal  coasting  into  maneuvering  targets  capable  of  tank- 
type  accelerations  are  quite  insignificant.  This  of  course,  is  based 
on  the  assumption  that  the  fissile  retains  its  last  correct  command 
prior  to  loss  of  track.  If  the  control  system  is  returned  to  its  neu¬ 
tral  positions,  the  miss  distances  increase  accordingly. 

Quantitative  data  on  the  missile  performance  requirements  to  over¬ 
come  several  important  perturbations  has  been  obtained.  A  sufficient 
variety  of  parameters  has  been  employed  in  an  attempt  to  cover  the  pro¬ 
bable  values  for  a  CCT-type  system.  It  has  been  shown  that  the  data 
obtained  can  be  extended  to  other  than  the  cases  shown,  provided  that 
•the  missile  velocity  and  engagement  range  are  such  that  tD  remains 
unchanged.  In  addition,  the  equations  and  analog  programs  for  each 
case  have  been  presented,  and  thus  further  data  could  be  obtained 
rapidly  if  necessary. 

Given  a  specific  set  of  missile  and  attack  parameters,  the  proper 
graphs  may  be  consulted  to  determine  the  maximum  acceleration  required 
to  overcome  given  disturbances.  For  example,  if  Vm  =  1000  ft/sec, 

Ro  =  6000  feet,  to/**"3 12  and  r  =  0.5.  What  missile  acceleration  is  re¬ 
quired  to  overcome  an  initial  leading  error  of  5°?  Graph  3  indicates 
that  50  ft/sec^  at  a  range  of  4800  feet  is  required. 
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Figure  I:  Definition  of  Angles 
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FIGURE  V:  Computer  Diagram  -  Forcing  Function  for  Case 
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TABLE  I 


Definitions  and  Symbols  for  Figure  II 


&1 »  ^2 •  k3 »  ^ 

Ao  flo  C0 

al»  bl 
ag,  b2 

1.  2.  3,  t 

%»  xt 

zm»  zt 

a 

♦ 

or 

Clo 


Initial  range 

Range  at  time  t  for  reference  collision  course 

Actual  range  at  time  t 

Auxiliary  lines,  parallel  to  R0 

Collision  triangle  for  constant  bearing  course 

Reference  missile  ana  target  positions  at  time  t 

Actual  (perturbed)  missile  &  target  positions  at 
time  t 

Auxiliary  lines  perpendicular  to  R0 

Longitudinal  missile  and  target  perturbations 

Traverse  missile  and  target  perturbations 

Missile  lead  angle 

Target  aspect  angle 

Line  of  sight  perturbation  angle 

Initial  LOS  angle  from  reference 
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TABLE  III:  MISS  DISTAKCES  FOR  CASE 
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